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Systems made by to heavy quarks play an important role in several high en- 
ergy experiments. The diversity, quantity and accuracy of the data currently 
being collected is impressive and includes: data on quarkonium formation 
from BES at BEPC, E835 at Fermilab, KEDR (upgraded) at VEPP-4M, 
(— i | and CLEO-III, CLEO-c; clean samples of charmonia produced in B-decays, 

in photon-photon fusion and in initial state radiation, from the B-meson 
factory experiments, BaBar at SLAC and Belle at KEK, including the un- 
expected observation of large amounts of associated (cc)(cc) production; 
the CDF and DO experiments at Fermilab measuring heavy quarkonia pro- 
duction from gluon-gluon fusion in pp annihilations at 2 TeV; the Selex 
experiment at Fermilab with the preliminary observation of possible candi- 
dates for doubly charmed baryons; ZEUS and HI, at DESY, studying char- 
monia production in photon-gluon fusion; PHENIX and STAR, at RHIC, 
and NA60, at CERN, studying charmonia production, and suppression, in 
heavy-ion collisions. In the near future, even larger data samples are ex- 
pected from the BES-III upgraded experiment, while the B factories and 
the Fermilab Tevatron will continue to supply valuable data for several 
years. Later on, new facilities will become operational (LHC at CERN, 
Panda at GSI, hopefully a Super-B factory, a Linear Collider, etc.) offering 
fantastic challenges and opportunities in this field. A comprehensive review 
of the experimental and theoretical status of heavy quarkonium physics 
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may be found in the Cern Yellow Report prepared by the Quarkonium 
Working Group 1 . See also the talks at the last QWG meeting at BNL (cf. 
http: / / www.qwg.to.infn.it /) . 

On the theory side, systems made by two heavy quarks are a rather 
unique laboratory. They are characterized by the existence of a hierarchy 
of energy scales in correspondence of which one can define a hierarchy of 
nonrelativistic effective field theries (EFT), each EFT has less degrees of 
freedom left dynamical and is simpler. Some of these physical scales are 
large and may be treated in perturbation theory. The occurrence of these 
two facts makes two heavy quark systems accessible in QCD. In particular 
the factorization of high and low energy scales realized in the EFTs allows 
us to study low energy QCD effects in a systematic and under control way. 
Ultimately the simplest EFT defined at the ultrasoft energy, pNRQCD, 
allows us to study the Quantum Mechanics of a non-Abelian field theory. 

2. Scales and EFTs 

The description of hadrons containing two heavy quarks is a rather chal- 
lenging problem, which adds to the complications of the bound state in 
field theory those coming from a nonperturbative QCD low-energy dynam- 
ics. A simplification is provided by the nonrelativistic nature suggested by 
the large mass of the heavy quarks and manifest in the spectrum pattern. 
Systems made by two heavy quarks are thus characterized by three energy 
scales, hierarchically ordered by the quark velocity v -C 1: the mass m 
(hard scale), the momentum transfer mv (soft scale), which is proportional 
to the inverse of the typical size of the system r, and the binding energy 
mv 2 (ultrasoft scale), which is proportional to the inverse of the typical 
time of the system. In bottomonium v 2 ~ 0.1, in charmonium v 2 ~ 0.3. In 
perturbation theory v <~ a s . Feynman diagrams will get contributions from 
all momentum regions associated with these scales. Since these momentum 
regions depend on a s , each Feynman diagram contributes to a given ob- 
servable with a series in a s and a non trivial counting. For energy scales 
close to Aqcd, the scale at which nonperturbative effects become dominant, 
perturbation theory breaks down and one has to rely on nonperturbative 
methods. Regardless of this, the non-relativistic hierarchy m ^> mv 3> mv 2 
will persist also below the Aqcd threshold. 

The wide span of energy scales involved makes also a lattice calculation 
in full QCD extremely challenging. We may, however, take advantage of 
the existence of a hierarchy of scales by substituting QCD with simpler but 
equivalent EFTs. A hierarchy of EFTs may be constructed by systemati- 
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cally integrating out modes associated to energy scales not relevant for the 
two quark system. Such integration is made in a matching procedure that 
enforces the equivalence between QCD and the EFT at a given order of 
the expansion in v and achieves a factorization between the high energy 
and the low energy contributions. By integrating out the hard modes one 
obtains Nonrelativistic QCD. In such EFT soft and ultrasoft scales are left 
dynamical and still their entanglement complicates calculation and power 
counting. We will focus here on the simplest EFT you can write down for 
two heavy quark systems, where only ultrasoft degrees of freedom remain 
dynamical. This is potential NRQCD a . 

3. pNRQCD 

pNRQCD 5-7 is the EFT for two heavy quark systems that follows from 
NRQCD by integrating out the soft scale. Here the role of the potential 
and the quantum mechanical nature of the problem are realized in the 
fact that the Schrodinger equation appears as zero problem for two quark 
states. We may distinguish two situations: 1) weakly coupled pNRQCD 
when mv 3> Aqcd, where the matching from NRQCD to pNRQCD may 
be performed in perturbation theory; 2) strongly coupled pNRQCD when 
mn <~ Aqcd, where the matching has to be nonperturbative. Recalling that 
r _1 ~ mv, these two situations correspond to systems with inverse typical 
radius smaller than or of the same order as Aqcd • 

3.1. Weakly coupled pNRQCD 

The effective degrees of freedom arc: low energy QQ states (that can be 
decomposed into a singlet and an octet field under colour transformations) 
with energy of order Aqcd, mv 2 and momentum p of order mv, plus ultra- 
soft gluons with energy and momentum of order Aqcd , mv 2 . All the gluon 
fields are multipole expanded (i.e. expanded in r). The Lagrangian is then 
given by terms of the type 



where the potential matching coefficients V n encode the non-analytic be- 
haviour in r. At leading order in the multipole expansion, the singlet sec- 
tor of the Lagrangian gives rise to equations of motion of the Schrodinger 



c k {m,u) 
m k 



x V n (ru',ru) x O n (u', mv 2 , A QCD ) r n . 



(1) 



a for an alternative and equivalent EFT (in the case in which Aqcd is the smallest scale) 
4 

see . 
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type. Each term in the pNRQCD Lagrangian has a definite power count- 
ing; retardation (or non-potential) effects start at the NLO in the multipolc 
expansion and are systematically encoded in the theory, they are typically 
related to nonperturbative effects 6 ' 7 . Resummation of large logs (typically 
logs of the ratio of energy and momentum scales) can be obtained using 
the renormalization group (RG) adapted to the case of correlated scales 
12 ; Poincare invariance is not lost, but shows up in some exact relations 
among the matching coefficients 9 . The renormalon subtraction may be 
implemented systematically. 7 
Applications of weakly coupled pNRQCD 

QCD Singlet Static potential. The singlet and octet potentials are well 
defined matching coefficients to be calculated in the perturbative matching. 
In 11 a determination of the singlet potential at three loops leading log has 
been obtained and correspondingly also a determination of ay showing how 
this quantity starts to depend on the infrared behaviour of the theory at 
three loops. The perturbative calculation of the static potential at (almost) 
three loops and with the RG improvement has been compared to the lattice 
calculation of the potential and found in good agreement up to about 0.25 
fm 17 . 

b and c masses. 

Heavy quarkonium is one of the most suitable system to extract a precise 
determination of the mass of the heavy quarks b and c. Perturbative de- 
terminations of the T(IS') and J/ip masses have been used to extract the 
b and c masses. The main uncertainty in these determinations comes from 
nonperturbative contributions (local and nonlocal condensates 6 ) together 
with possible effects due to subleading renormalons 1 . A recent analysis 
performed by the QWG 1 and based on all the previous determinations 
indicates an error of about 50 MeV both for the bottom (1% error) and in 
the charm (4% error) mass. 
Perturbative quarkonium spectrum. 

B c mass. Table 1 shows some recent determinations of the B c mass in 
perturbation theory at NNLO accuracy compared with a recent lattice 
study 16 and the value of the CDF experimental B c mass. This would 
support the assumption that nonperturbative contributions to the quarko- 
nium ground state are of the same magnitude as NNLO or even NNNLO 
corrections, which would be consistent with a mv 2 > Aqcd power count- 
ing. Hyperflne splittings, cc, bb, B c ground state hyperfine splittings 
have been recently calculated at NLL in 19 . The prediction for r]b mass is 
M(r)b) = 9421 ± 10 (th) ± 9 (Sa s ) MeV. The logs resummation seems to be 
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important. If the experimental error in future measurements of M(rjb) will 
not exceed few Mev, the bottomonium hyperfine separation will become a 
competitive source of a s (Mz) with an estimated accuracy of ±0.003. 
Radiative transitions (Ml). A theory of Ml transitions in heavy quarko- 
nium has been recently formulated using pNRQCD 20 . This may shed some 
light on recent CLEO results on radiative Ml transitions in the rjb search 
that have ruled out several models. No large anomalous quarkonium mag- 
netic moment is generated. 

Seminclusive radiative decays of T(15). In 22 the end-point region of 
the photon spectrum in semi-inclusive radiative decays of heavy quarko- 
nium has been discussed using Soft-Collinear Effective Theory and pN- 
RQCD. Including the octet contributions a good understanding of the ex- 
perimental data is obtained. 

Gluelump spectrum. In pNRQCD 6 ' 21 the full structure of the gluelump 
spectrum has been studied, obtaining model independent predictions on 
the shape, the pattern, the degeneracy and the multiplct structure of the 
hybrid static energies for small QQ distances that well match and interpret 
the existing lattice data. 

Properties of baryons made of two or three heavy quarks. Recently 
the SELEX experiment has detected first signals from three-body bound 
states made of two heavy quarks and a light one. The two heavy quark part 
of such systems may be treated in pNRQCD and a precise predication for 
the hyperfine interaction may be obtained 23 . 



B c mass (McV) 


18 (expt) 


16 (lattice) 15 (NNLO) 13 (NNLO) 14 (NNLO) 


6287 ± 4.8 ± 1.1 


6304 ±12^ 2 6326(29) 6324(22) 6307(17) 



3.2. Strongly coupled pNRQCD 



In this case the matching to pNRQCD is nonperturbative 25 . In the situa- 
tion where the other degrees of freedom (like those associated with heavy- 
light meson pair threshold production and heavy hybrids) develop a mass 
gap of order Aqcei, the quarkonium singlet field S remains as the only low 
energy dynamical degree of freedom in the pNRQCD Lagrangian (if no 
ultrasoft pions are considered), which reads 25 > 26 ' 6 > 7 : 

ipNRQCD - Tr {St (ido - |^ - Vs(r)) s}. (2) 

In this regime we recover the quark potential singlet model from pNRQCD. 
The matching potential Vs (static and relativistic corrections) is nonper- 
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turbative: the real part controls the spectrum and the imaginary part con- 
trols the inclusive decays. The potential is calculated in the nonperturba- 
tive matching procedure between NRQCD and pNRQCD 25 ' 7 . Advantages 
of this approach include: factorization of hard (in the NRQCD matching 
coefficients) and soft scales (contained in Wilson loops or nonlocal gluon 
correlators); the low energy objects being only glue dependent, confine- 
ment investigations, on the lattice and in QCD vacuum models become 
feasible 24 ; the existence of a power counting indicating leading and sub- 
leading terms in quantum-mechanical perturbation theory; the quantum 
mechanical divergences (like the ones coming from iterations of spin delta 
potentials) are absorbed by NRQCD matching coefficients. The potentials 
evaluated on the lattice once used in the Schrodingcr equation produce the 
spectrum. The calculations involve only QCD parameters (at some scale 
and in some scheme). 

Applications of strongly coupled pNRQCD 

Nonperturbative potentials and Spectrum. Recently the multilevel 
algorithm has been applied to the lattice evaluation of field strength inser- 
tion inside the Wilson loop average, producing very precise data for the spin 
dependent potentials and a first evaluation of the nonperturbative potential 
at order 1/m 27 . This is the first step towards a precise determination of 
the nonperturbative matching potentials on the lattice. 
Decays. The inclusive quarkonium decay widths in pNRQCD can be 
factorized with respect to the wave function (or its derivatives) calcu- 
lated in zero, which is suggestive of the early potential models results: 
T(H — > LH) = F(a s , Aqcd) • IV^O)! 2 - Similar expressions hold for the elec- 
tromagnetic decays. However, the coefficient F depends here both on a s 
and Aqcd- In particular all NRQCD matrix elements, including the octet 
ones, can be expressed through pNRQCD as products of universal non- 
perturbative factors by the squares of the quarkonium wave functions (or 
derivatives of it) at the origin. The nonperturbative factors are typically 
integral of nonlocal electric or magnetic correlators and thus depends on 
the glue but not on the quarkonium state 26 . Typically F contains both the 
NRQCD matching coefficients at the hard scale m and the nonperturbative 
correlators at the low energy scale Aqcd- The nonperturbative correlators, 
being state independent, are in a smaller number than the nonperturba- 
tive NRQCD matrix elements and thus the predictive power is increased 
in going from NRQCD to pNRQCD. Thus, having fixed the nonpertur- 
bative parameters on charmonium decays, new model-independent QCD 
predictions can be obtained for the bottomonium decay widths 26 . 
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